Exact and approximate distributions for the product of Dirichlet components

It is well known that X/(X -f Y) has the beta distribution when X and Y follow the Dirichlet distribution. Linear combinations of the form aX + pY have also been studied in Provost and Cheong [24]. In this paper, we derive the exact distribution of the product P = XY (involving the Gauss hypergeometric function) and the corresponding moment properties. We also propose an approximation and show evidence to prove its robustness. This approximation will be useful especially to the practitioners of the Dirichlet distribution.